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Abstract 

In relativistic kinetic theory, which underlies relativistic hydrodynamics, the molecular chaos 
hypothesis stands at the basis of the equilibrium Maxwell-Jiittner probability distribution for 
the four- momentum p a . We investigate the possibility of validating this hypothesis by means 
of microscopic relativistic dynamics. We do this by introducing a model of relativistic colliding 
particles, and studying its dynamics. We verify the validity of the molecular chaos hypothesis, and 
of the Maxwell-Jiittner distributions for our model. Two linear relations between temperature and 
average kinetic energy are obtained in classical and ultrarelativistic regimes. 
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I. INTRODUCTION 



The study of relativistic fluids^both from the hydrodynamic and kinetic point of view has 
been widely investigated . In this context, the relativistic Boltzmann equation 



„ df dfF» r fl , dy d 3 P * dy, h s 



represents the best known tool, which is based on a molecular chaos hypothesis, like the 
Boltzmann equation in classical kinetic theory. Here x v , p u , F v are respectively the position, 
momentum and force four-vectors, m is the rest mass, Q is the interaction cross-section, 
and / is the single particle distribution function. Collisionless relativistic plasmas are inves- 
tigated by means of the relativistic Vlasov equation, obtained neglecting the collision term 
in Eq. ([T|) . The equations of relativistic Hydrodynamics, which macroscopically describe rel- 
ativistic fluids, are derived also from Eq.l f T j) . similarly to the classical case. 
The chaotic hypothesis, which underlies Eq. fTjl . explains how the microscopic components 
of a fluid reach a local equilibrium state. Classically, it is well established that this is a con- 
sequence of the interactions among the particles, as illustrated, for instance, by molecular 
dynamics 

In order to investigate the validity of the molecular chaos assumption in relativistic kinetic 
theory, we propose a simple model of N relativistic colliding particles, and investigate the 
properties of its dynamics. 

In fact, to the best of our knowledge, many particle relativistic systems have only been stud- 
ied either from a kinetic or hydrodynamic point of view, because the microscopic dynamics 
of such particle systems presents many difficulties. For instance, it is highly problematic 
to write covariant hamiltonians (and the related 4- vector equations of motion) for the sys- 
tems. Other difficulties concern: the choice of the reference frame, since every particle has 
a different proper time; the form of the interaction potential, since the action and reaction 
principle holds only for contact interactions; the effects of length contraction and time dila- 
tion. The consequence of this is that, as far as we know, no direct microscopic evidence for 
the molecular chaos hypothesis in relativistic dynamics has been provided. 
To overcome this difficulty, we propose a non-covariant hamiltonian written with respect to 
the center of mass frame, taken as the Lorentz rest frame, which yields the non-covariant 
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equations of motion 



dxj 



dt 
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where N is the number of particles. For the force F^ca wg p r0 p 0se ^ Q use 
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where = (r^ — rj), r^j = \rij\, ^^ CA is the Weeks- Chandler- Andersen interaction poten- 
tial [7]; the quantities e and a are obtained from the Lennard- Jones (LJ) potential which 
defines ^^ CA , and represent respectively the depth of the LJ potential, and the distance at 
which it changes sign. 

Therefore, particles move according to the relativistic dynamics when they do not interact, 
while their interactions are modelled classically, so that the total momentum and the total 
kinetic energy of particles are preserved by the collision process. Although this is not com- 
pletely rigorous, our procedure meets all the microscopic requirements of relativistic kinetic 
theory, i.e. the invariance of the momentum 4-vectors. 

In this paper we simulate a 2D system of A^ relativistic particles (with A^ = 28), through 
a MD algorithm, which implements the equations of motion (|2I3|) with periodic boundary 
conditions, for a density p = N/A = 0.2 (with A the cell area), which is not a low density 
case. The simulations are performed for different initial kinetic energies corresponding to 
classical, relativistic and ultrarelativistic regimes. Furthermore, we take e = a = 1. 
In the low density limit, the contribution of the collisions is expected to become negligible, 
and the dynamics to tend to a fully covariant dynamics. 



II. RESULTS 



Our results show that the simulated systems all reach an equilibrium state since their 
observables, such as the pressure, converge to an equilibrium value, while the probability 
distribution functions (PDFs) of the values of microscopic quantities like momentum p x and 
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kinetic energy £ reach an invariant form. In particular, we find that the PDFs of p x reduce 
to the Maxwell-Boltzmann (MB) distribution in the classical limit, as desired. This is due 
to chaos in the dynamics, which is evidenced by the fact that the numerically evaluated 
largest Lyapounov exponents are positive. 



A. Probability Distribution Functions 

The standard relativistic kinetic theory predicts that the PDF of p a has the form of 
the Maxwell- Jiittner (MJ) distribution, f MJ = dexp (— [7 Q p a //c B T),with d a normalization 
constant and U a the hydrodynamic four- velocity (with U z = 0) Q, Q|. In the local rest 
frame, fuj can be written as 



fMj(px,P y ) = d exp y~ a ]J + X j ' ( 4 ) 

where p x , p y are the spatial components of p a , c is the speed of light, and where d and a are 
two constants related by the normalization condition 

d= G^p) TTZ • (5) 

As well known , a involves the temperature of the system [llj], because 

a = m c 2 /ksT . (6) 
Integrating Eq.(0J) over p y , one obtains the PDF for p x only: 



g MJ (p x )=2m cd x l-% 5 + l -KxlaJ-^ + l) , (7) 

where Ki(x) is the modified K-Bessel function of first order. Considering the kinetic energy 
£ = ca/p 2 + mpC 2 — m c 2 , Eq.(j3J can also be rewritten as 

h M j{i) = ^rd (£ + m c 2 ) • exp (-a + m ° C j . (8) 
<r \ m c^ / 

It is interesting to observe that, if an expression like the MB distribution was written for 
the relativistic p x , i.e. if one started from 

fMB(Px,Py) = fMB(Px)fMB(Py) = exp ( -O — T ^ ) , (9) 



the PDF of the relativistic kinetic energy £, after some calculations, would take the form 



hMEs(0 



2ae a 
mlc 4 



(£ + m c 2 ) exp 



~ (£ + ™oC : 
-a— 



,2\2 



mgC 4 



(10) 



Comparing Eq.s (|91 10|) with Eq.s (|7I 8|) . one notices that the MJ distribution is not merely 
the MB distribution with the relativistic p x and £ in place of the classical momentum and 
kinetic energy. 

We fit the histograms constructed through our MD simulations to the PDFs given above, 
and for simplicity we take itlq — c — 1. 

The following figures are obtained for different mean kinetic energies, where the mean kinetic 
energy is the time average of the total kinetic energy divided by the number of particles. 
The histograms are constructed recording the instantaneous values of momentum p x and 
kinetic energy £ for a given particle. This operation is repeated every 200 timesteps, in 
order to decorrelate the recorded data. 



Mean kinetic energy per particle = 9.87 x 10 
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FIG. 1: Fit of the data for momentum p x on a log scale (left panel). Fitted histograms of the kinetic energy f in linear scale 
(right panel). The parameter of the Maxwell- Jiittner PDFs takes the value a = 10.4875 for the PDF ofp^ yielding kgT = 0.095 
(left panel), and a = 11.2595 for the PDF of £ leading to k B T = 0.089 (right panel). The classical Maxwell Boltzmann PDFs 
fits well the data only in the low energy cases. 
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Mean kinetic energy per particle = 9.83 x 10 
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FIG. 2: Fit of the data for momentum p x on a log scale (left panel). Fitted histograms of the kinetic energy £ in linear scale 
(right panel). The parameter of the Maxwell-Jiittner PDFs takes the value a = 1.4067 for the PDF of p x yielding fcsT = 0.711 
(left panel), and a = 1.4225 for the PDF of £ leading to k B T = 0.703 (right panel). The classical MJ distributions fits better 
the data than the MB ones at these energies. 




FIG. 3: Fit of the data for momentum p x on a log scale (left panel). Fitted histograms of the kinetic energy £ in linear scale 
(right panel). The parameter of the Maxwell-Jiittner PDFs takes the value a = 0.2540 for the PDF of p x yielding k B T = 3.937 
(left panel), and a = 0.2501 for the PDF of £ leading to k B T = 3.998 (right panel). The MB docs not fit the data at these 
ultrarclativistic energies. 



6 



For the Maxwell- Jiittner PDFs, if the parameters a and d are obtained as independent 
parameters by fitting the numerical data to Eq.s (|7|8|) . the normalization condition (jSJ) both 
for p x and £ is verified. This indicates that the MJ-PDF is indeed appropriate for our data, 
and that the data are consistent. 



B. Measurement of temperature for a relativistic system 



The microscopic definition of the temperature of a system composed by relativistic parti- 
cles is an open issue o| . However, the Maxwell- Jiittner PDF contains one parameter, which, 
in analogy with the classical Maxwell-Boltzmann PDF is identified with the quantity /cbT. 
Therefore, observing that the Maxwell- Jiittner PDFs fit well our histograms, it becomes 
reasonable to assume for our system a -1 = k B T as a definition of temperature obtained 
from the microscopic dynamics. 
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FIG. 4: Plot of temperature vs mean kinetic energy per particle. Temperatures are calculated through Eq. from the 
fitting parameter a -1 = kgT, obtained both for the kinetic energy(*) and the p x (o). The two different calculations yield 
indistinguishable values in this figure. In the left panel the classical regime is plotted, while in the right one the ultrarclativistic 
limit is shown. In both panels the dotted lines represent the relation between kgT and £ valid in the low energy cases; the 
dash-dotted lines the relation between fcgT and § valid in the high energy cases. 



A linear relation between this temperature and the mean value of the kinetic energy per 
particle has been found for the classical and ultrarelativistic cases. For k B T = a^ 1 < 0.1 
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(classical regime), the relation was found to be, as expected, ksT = £, while for ksT = 
a^ 1 > 1 (relativistic and ultrarelativistic regimes), we verified a linear relation of the form 
/cgT = 0.56£ + 0.21. The transition between the two regimes takes place in a small range of 
kinetic energy values. 

III. CONCLUSIONS 

In this paper we have tested a 2D molecular dynamics model intended to simulate the 
microscopic dynamics of N relativistic colliding particles, with total constant energy E, and 
have observed its relaxation to an equilibrium state. 

Our model satisfies the requirements of momentum and kinetic energy conservation before 
and after the collisions, underlying the equilibrium relativistic kinetic theory. The histograms 
found by these simulations for the momentum p x , and for the kinetic energy £ are well fitted 
by the PDFs of the standard relativistic kinetic theory, i.e. by the PDFs derived from the 
MJ distributions. In addition to this, the statistics of the dynamics of our model reduces to 
the classical one when the kinetic energy takes small values. 

Our model suffers from the difficulties of not being fully relativistic, because the particle 
interactions are treated classically; therefore, it becomes more and more acceptable as the 
particle density decreases, or the collision rate tends to zero making the dynamics tend 
to a fully covariant form. Moreover, as we are going to report in (ill ], reducing densities 
does not produce any qualitatively different result, which indicates that in the limit of 
low collision rates the macroscopic behaviour of our systems is not substantially different 
from that of the higher density cases. This, together with the observed validity of the MJ 
distributions, provides a justification for our model, as a tool to simulate relativistic many 
particle systems. Otherwise, if this model is accepted, it affords a microscopic justification 
of the relativistic molecular chaos hypothesis, underlying relativistic kinetic theory and 
relativistic hydrodynamics. 

Furthermore, linear relations of temperature and mean kinetic energy have been found both 
in classical and ultrarelativistic regimes. This allows us to obtain a definition of temperature 
in a relativistic system, something rather problematic in general Q] , which deserves further 
investigations. 
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